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Abstract 

We compute moving eigenvalues and the eigenvectors of the linear re- 
normalization group transformation for observables along the renormalized 
trajectory of the hierarchical non-linear 0(iV)-invariant cr-model by means of 
perturbation theory in the running coupling constant. Moving eigenvectors 
are defined as solutions to a Callan-Symanzik type equation. 



1 Introduction 



Hasenfratz and Niedermayer calculated in ||HN|| the asymptotic form of the renor- 
malized trajectory called perfect action of the two dimensional 0(N)-model with 
non-linear block spin. They observed that the block spin transformation could be 
treated classically in the asymptotically free domain to an excellent approximation. 
In ||WX|| we introduced an improvement scheme in which the zeroth order approxi- 



mation is the perfect action. This scheme is a combination of perturbation theory in 
the running coupling with the idea of scaling. We worked it out for the 0(iV)-model 
with linear block spin in the hierarchical approximation. The improved action is 
not yet the complete improved theory. The complete improved theory should also 
include improved observables. In this paper we show how the improvement program 
in [WX] can be extended to observables for the hierarchical 0(iV)-model. 

In Wilson's renormalization group [|WK|| a Euclidean quantum field theory is 
thought of in terms of the flow of effective interactions V(4>) generated by a block 
spin transformation TZ. In the following <f> will denote an iV-component real valued 
field on a two dimensional unit lattice. The interactions will be invariant under 
global O(N) transformations. TZ will respect global O(N) invariance. The fixed 
points IZV* — K of TZ describe scale invariant (continuum) field theories. The 
eigenvalues of the linearized block spin transformation at a fixed point determine 
the associated critical exponents. The two dimensional 0(iV)-invariant non-linear 
cr-model is expected to possess only two fixed points for N > 2 (and N > 1 in the 
hierarchical approximation), an asymptotic, unstable ultraviolet fixed point Vuv 
and a stable infrared fixed point Vir. The model is asymptotically free in the 
ultraviolet: Vuv is a theory of iV — 1 massless free fields. In the hierarchical case 
Vuv is better thought of from the orthogonal point of view as a high temperature 
fixed point of a single component radial field. The fixed points are expected to 
be connected by a one dimensional renormalized trajectory which is stable under 
block spin transformations. The points on the renormalized trajectory describe 
(continuum) theories with broken scale invariance. Interactions on the renormalized 
trajectory are said to scale. The renormalized trajectory is both of principal and 
of practical interest since it admits the description of a true continuum theory (not 
a lattice approximation) in terms of a lattice theory. It is the ultimate improved 
action in the sense of Symanzik JS|. 

Let the renormalized trajectory be parametrized in terms of a local coordinate 
/. Recall that it is a curve in the space of O(N) invariant interactions. The block 
spin transformation then becomes 

KVktW) = VRrWif))- (1) 
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A natural local coordinate in the vicinity of the asymptotically free fixed point is the 
running coupling. The running coupling will be defined as the leading trilinear o- 
7r-7r-vertex. In the hierarchical approximation it is identical with the inverse radius. 
The complete dynamics of the renormalization group on the renormalized trajectory 
is encoded in the flow of the running coupling. We call /3(f) Callan-Symanzik 
function since it defines this discrete flow. 

A block spin transformation comes together with a linear mapping of observables 
which preserves their expectation values. This mapping is the linearized renormali- 
zation group transformation. We shall consider this linear mapping CTZ to the block 
spin transformation over the renormalized trajectory. An observable will be called 
an eigenvector with eigenvalue e(/) if it satisfies the equation 

c RT no(<p\f) = <P(f))o(mf))- (2) 

In accordance with the terminology for the interaction we will speak of this property 
as scaling. The space of eigenvectors defines a moving frame in the tangent space 
over the renormalized trajectory. The set of eigenvalues will be called the spectrum 
of the theory. It is exactly calculable at the fixed point. Moving away from the fixed 
point on the renormalized trajectory it becomes perturbed. In the vicinity of the 
ultraviolet fixed point it is therefore natural to perform a perturbation expansion 
in the running coupling to determine the spectrum from (^]). In this paper we 
address the question of computing these perturbative corrections. At the fixed 
point (Ei) becomes the eigenvalue equation for the critical exponents. Away from 
the fixed point the spectrum governs the rate of attraction towards the renormalized 
trajectory. The set of eigenvectors forms a basis of the linear space of observables. 
Associated with this basis is a set of fusion rules. We will give a general expression 
for correlation functions of eigenvectors on the renormalized trajectory in terms of 
the spectrum and the fusion rules. The investigation of the spectrum and the moving 



frame complete the analysis begun in [ |WX|| . 



We should mention that the two dimensional hierarchical 0(iV)-model has been 
rigorously studied by Gawedzki and Kupiainen in |GK|, where both the existence 



of an ultraviolet limit and the existence of a renormalized trajectory has been es- 
tablished beyond perturbation theory. This result has also been obtained by Pordt 
and Reisz in ||PR[1 using a slightly different technology. Both groups however do not 
address the question of the spectrum. A rigorous study of the full 0(N)-mode\ both 
in the ultraviolet and more importantly in the infrared domain is still an outstanding 
problem. 
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2 Hierarchical 0(N) model 



Let us consider the renormalization group flow generated by the hierarchical block 
spin transformation^ 



e 



in D = 2 dimensions with scale parameter L = v2. The interaction V((f>) is taken 
to be a function of a real ^-component variable (p. 

d^ 7 (C) = (2 7 r 7 )^e^d 7V C (4) 

is the Gaussian measure on 1R N with mean zero and covariance diag(7, . . . , 7). The 
subspace of O(N) invariant interactions is stable under (|3|). We will restrict our 
attention to this subspace. O(N) invariance requires V(<p) to be a function of the 
modulus \<p\. The normalization constant in (|3]) is conveniently chosen such that the 
interaction is always zero at its minimum. We will denote this minimum by r and 
call it the radius of the potential. 

Let us define the renormalized trajectory as the curve Vrt(^|/) in the space 
of O(N) invariant potentials, parametrized by the inverse radius / = -, with the 
following two properties: 

1) Vrt(0|/) is stable under the block spin transformation 71. It follows that there 
exists a function (3(f) such that 

KV RT (cf>\f) = V RT (<P\(3(f))- (5) 

In other words a block spin transformation acts on the renormalized trajectory by 
a transformation of the coordinate given by a /3-function. 

2) The asymptotic behavior of Vrt(^|/) as the running coupling / goes to zero is 
given by 



VrtW) = vj£(<p\f) + o(f 



1 2 1 

2 



v$W) = ^(M-yJ • ( 6 ) 

Up to corrections of second order in the running coupling the renormalized trajectory 
coincides with the perfect action. 



1 Note that (|^) differs from the transformation used in [ WX ] . We have exchanged the integration 
with the rescaling step. 
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At least in terms of perturbation theory in the running coupling (|5|) and (||) have a 
unique solution as has been shown in |[WX|| . Let us remark that (|[) is compatible 
with (|) since it reproduces itself under (|3|) according to 

KVj&Wf) = V$(mf)) + 0(f), 

(3(f) = f + 0(f). (7) 

up to second order corrections. The general form of the solution to (^) and (|5[) in 
terms of perturbation theory is 



n=2 
1 



Pz{f) — n 1" ^ C 2,2m+2/ 2m+2 : 

^7 m =0 



oo 

p2n+m 



Pn{f) — X/ C n,2m+nf^ 
m=0 

oo 



m=0 

The coefficients up to fifth order in the running coupling turn out to be given by 

3(N- 1) 



C2,2 



4 ' 

'9(iV-l) 2 61(JV-1) 



C2 ' 4 = 8 28— 7 ' 



C3,3 



7(JV- 1) 



18 

5(iV-l) 2 257(iV-l) 



C3 ' 5 = 8 + 180 7 ' 



C4,4 
C5,5 
&3 



15(iV 



56 

31(iV- 1) 

150 
(iV-l)7 

2 ' 



6s = /^ + ^v_i)A 2) 



(9) 
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It follows that the model is asymptotically free for N > 1. It is custom to use g = f 2 
instead of / as running coupling in discussions on asymptotic freedom. From the 
/3-function to fifth order in / it follows that the recursion relation to third order in 
g is given by 

oo 
m=2 

b 2 = (iV -1)7, 

Thus g is marginally relevant for N > 1. 

This definition of a renormalized trajectory does not refer to a continuum limit 
procedure. It is nevertheless identical with the continuum limit effective potential 
of models in the O(N) universality class in the hierarchical renormalization scheme. 
To prove this we can perform the continuum limit using @ as bare interaction. 
Define the bare coupling to be the n-fold preimage (3~ n (f) of a renormalized value 
/. By construction of the renormalized trajectory it follows that 

n n V RT {^- n {f)) = Vartflf) (11) 

for all numbers n of renormalization group steps. The continuum limit n — > oo is 
immediately performed since the right hand side of (|TT|) is independent of n. All 
that is needed is an analysis of the recursion relation defined by the /3-function or 
rather its inverse. This is a comparatively easy task. See for instance ||GK]|. For 



AT > 1 it follows from fllPf) that the bare coupling tends to zero as n goes to infinity. 
(One needs to take into account logarithmic corrections piled up by the term of 
third order in g, fifth order in /.) It follows that the perturbation expansion in 
the running coupling is valid. To make contact with the hierarchical real world one 
should also assign a scale, for instance in form of a lattice spacing a, to the point on 
the renormalized trajectory where the running coupling is given by the renormalized 
value /. The bare cutoff in (|TT|) is then L~ n a with L the block scale. 

The sceptical reader may worry to what extent this construction is connected 
with the continuum limit of a bare theory defined by his favorite O(N) invariant 
interaction. Consider for instance the standard interaction of the linear 0(N)-model 
defined by 

V{<P\f) = V 2 (f-j)j • (12) 
Its continuum limit is constructed as the result of the infinite iteration 

lim TZ n V(^ n (f)) = v cor Mf)- (is) 



5 



Here the bare coupling f- n (f) is tuned such that the minimum of H n V((j)\f- n (f)) 
is located at the renomalized radius \(f>\ = j. The continuum limit is universal as we 
have learnt from the work of Wilson ||WK |. Therefore the connection to the above 



definition is simply 

V cont (<P\f) = V RT W)- (14) 
A rigorous proof of the existence of and, inbetween the lines, also of ( |H|) has 
been given by Gawedzki and Kupiainen ||GK |. More generally, the set of all bare 



interactions sharing (|8|) as their common continuum limit defines the universality 
class of the hierarchical 0(iV)-model. (|12|) is known to belong to this class and so 
is the original model with sharp constraint. So what we do, when we define the 
renormalized trajectory by the above two conditions, is to completely disentangle 
the admittingly also interesting question if, how, and at what pace the continuum 
limit is reached by some particular bare model. 

In practice we may not be able to compute (f|) to all orders of perturbation theory 
in the running coupling for the models we are really interested in, for example the 
full nonlinear 0(iV)-model in terms of block spin transformation on a unit lattice. 
In this situation the best one can do is to take the highest order approximation to 
([§D accessible as bare interaction. In the hierarchical model for instance already the 
second approximant 

turns out to be an excellent starting point for a numerical study of the renormalized 
trajectory. Its main property is of course 

= V$(mf)) + 0(P(f) 3 )- (16) 

It therefore coincides with the renormalized trajectory up to corrections of third 
order in /. Let us finally mention that the third order approximant would be ideally 
suited for a rigorous construction along the lines of Gawedzki and Kupiainen 



and of Pordt and Reisz ||PR]| . The perturbative part is trivial since the action 
reproduces itself to third order. The construction therefore reduces to the proof of a 
stability bound which controls non perturbative corrections. This bound is already 
implicitly contained in the rigorous work on the 0(iV)-model with standard bare 
interaction (|12|). 



3 Perfect Observables 

Let us consider the hierarchical renormalization group transformation for local ob- 
servables corresponding to (H) in the case of D = 2 dimensions with scale parameter 
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L = y/2. It is given by the linear transformation 

ZvWW ~ /(WC )e-^+C) • 

The block volume is L D = 2. ( |TTD is the linearization of (H) divided by the block 
volume 

CvTWty) = \%-n{V + zO)\ z=0 . (18) 
2 az 

In other words the transformation of observables is the tangent map of the trans- 
formation of the potential. 

The setup of perturbation theory for the linear transformation fllTD is as follows. 
Since the model is 0(iV)-invariant we can choose the block spin to be given by 
ifi = (r + \l/)e with e any N component unit vector. That is, we can trade for 
\& = l^l — r. Observables will also be assumed 0(iV)-invariant. The shift by r serves 
to place us into the minimum of the potential. We then decompose orthogonally 
the fluctuation field into a radial and a tangential part with respect to the direction 
of the block spin ip. The decomposition is £ = ae + n. The one component variable 
a is the radial fluctuation field. The (N — l)-component variable 7r is the tangential 
fluctuation field. The Gaussian measure factorizes into 

dfjL y (Q = d/i 7 (cr) d / u 7 (7r). (19) 

We then insert this decomposition of £ in the formula for the potential on the 
renormalized trajectory and extract the terms of order zero in the running coupling. 
The result is 

V RT ((r + V + a)e + ir\f) = (*±f£ + y$(¥ + a, ir\f). (20) 

27 

The zeroth order terms cannot be treated as perturbations. They are fortunately 
only linear and quadratic in the radial fluctuation field. They and can therefore be 
accomodated for using the identity 

-(w+<t) 2 1 _*£ 

d/i 7 (a)e *y =dyUi ( a + _)_ e 4 7 . ( 2 l) 

We then shift the radial fluctuation field to £ = a + -I. After these standard 
manipulations the transformation fllTD becomes 



£ RT fiO((r + vl/)e)= ^ yq>J ^ ^-T^T^ 22 
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At this point perturbation theory is applicable. The expansion parameter is the 
inverse radius f — -■ For the observables we have in mind we also have to perform 
an expansion in terms of the running coupling. 

An observable is called a moving eigenvector to the moving eigenvalue e(f ) if it 
satisfies the renormalization group equation 

Cnr-RO^f) = e(/3(/))O(0|/3(/)). (23) 

Solutions to (|23"D are perfect observables. A parametrization well suited for pertur- 
bation theory is 

OO / 1 \ n 

0W) = £ (M - 7) • (24) 

Here the coefficients Q n {f) are taken to be power series in the running coupling /. 
(The coefficients are not expected to exhibit any singularity at zero coupling.) We 
will organize the solutions of (|23|) according to their zeroth order coefficients. The 
zeroth order observables are simply normal ordered monomials. Let us also perform 
the perturbation expansion for the observables of the form (^4|). Define 

WW) = |0| - j. (25) 

As in the potential we find a term of order zero in the running coupling. Separating 
it off we obtain 

W((j + | + + *\f) = | + £ + ^ (1) (| + & tt|/). (26) 

In the limit when the radius becomes infinite the transformation flUTD reduces to 
the convolution with the radial Gaussian measure. This transformation is identical 
with that of a one component scalar field in two dimensions at the high temperature 
fixed point. We are therefore immediately lead to normal ordered monomials. Put 
j' — Then it follows that 

/^(0:(§+0 B V=^:*"V. (27) 

Therefore the eigenvectors are normal ordered monomials to zeroth order as ex- 
pected. The zeroth order spectrum is e n = Let us write the normal ordered 
monomials in the form 

n 

: V 1 : y = £ h n ^ m . (28) 

m=0 
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They are of course rescaled Hermite polynomials in the variable We then write 
the associated sequence of observables defined as solutions of the scaling equation 
in the form 



in 



O n {4>\f)= Y.Qn,m{f) (H-t) • (29) 
m=0 V J J 

The coefficients Q n ,m{f) are given by the power series expansions 

oo 

Qn, m (f)=E d S m f l (3°) 
1=0 

with zero order coefficients of the normal ordered form 

Qn]m = h"n,m- (31) 

That is, the observables are perturbations of normal ordered monomials. The per- 
turbative form of the moving spectrum is 

oo 

*»(/) = E (32) 

m=0 

with zero order coefficients 

€> = ±. (33) 

This completes the setup for the improvement program for observables. Let us then 
turn to the question of how to compute the higher coefficients in the expansions (|30| ) 
and fl32~|) for the moving eigenvectors and observables. The strategy is an adaptation 
of the improvement program for the potential. Let us choose the first (nontrivial) 
observable Oi((j)\f) as an example and perform a third order computation to some 
detail. The zeroth order argument provides us with the information that 

OMf) = WW) + O(f), (34) 
ex(f) = l + 0(f). (35) 

Let us denote the zeroth order approximation by 0(<f)\f) = W(<p\f). As we will 
see this observable already scales to second order. We immediately perform a per- 
turbation expansion to third order for the effective observable of the zeroth order 
approximation. The result is 

(£ RT n)ow) = <f){Qo{n+w{ct>\f)+Q2{r)w{ ( t ) \r) 2 )+o{{fr), 
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5{N- 1)7 



2 



(•/\3 



Qo(/0 = - v 3 y ' (/ 

W) = - (jV ~ 1)7 (/') 3 , (36) 

where the old running coupling / is expressed in terms of the new running coupling 
/' = 0(f)- (^rtT^- is understood as an operator which is applied to the function 
0(4>\f).) This change of coupling prepares in particular the ground for further 
iterations. Let us emphasize that not only the coefficients but also the coordinate 
functions depend on the running coupling and have to be adjusted. From this 
we conclude that the observable already scales to second order. The eigenvalue is 
therefore already correct to second order. To third order both a constant and a term 
quadratic in W((p\f) are generated. Therefore the observable does not reproduce 
its dependence on the field to third order. To find the observable which scales to 
third order we make the ansatz 

0(<j>\f) = d f + W(<j>\f) + d 2 fW{4>\ff- (37) 

The ansatz involves two improvement parameters do and d 2 . To determine their 
value one again computes the effective observable starting from (|37|) . The expansion 
gives 

(£stR) 0(01/) = £(/') (<?„(/') + WW) + <M/')»'M/') 2 ) + o((/') 4 ), 
Qo(f') = (- 5(iV 3 1)1-2 + dr, + 2d„) (/?, 

W) = (-^1 + I ) (/<)*. (38) 

The improvement parameters are then determined by the simple linear equations 

-(TV -1) 7 d 2 

d2 ~ — — + y 

-5(N - \W 

do = 1 3 )J +d 2l + 2d . (39) 



The solution is 



d 2 = (iV 1)7 , do = 2(N - 1) 7 2 . (40) 
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To third order we therefore find 

OMf) = 2(iV - 1) 7 2 / 3 + W(<f>\f) - (N f 3 W(<f>\f) 2 + 0(/ 4 ) (41) 
together with the eigenvalue 

ei(f) = \+ i -^ 1 f 2 + 0(n. (42) 

Note that the eigenvalue does not have a term of third order in the running coupling. 
This scheme is iterated in the obvious manner. Suppose that we have computed the 
observable to order s in the running coupling 

OMf) = 0[ s \4>\f) + 0{f s+l ). (43) 

Suppose that the order s improved observable is given by the general form 

oi s) W) = EQil(f)wWT, 

n=0 

Q$(f) = h 

Qii(f) = f 4lr. (44) 

r=n+l 

Then by induction it follows that the effective observable to order s + 1 is again of 
this general form. It follows that to every order of perturbation theory only finitely 
many powers of the coordinate functions W(<f>\f) appear. The sums turn out to 
involve either even or odd powers in the running coupling respectively. We assume 
that 0[ s \(J)\f) scales to order s. That is, 

(C RT TZ) 0[ s \<P\f) = 0{ s \4>\f) + 0((f ) s+1 ). (45) 
Then we take an ansatz of the form to order s + 1 treating the coefficients of 



order s + 1 in the running coupling as improvement parameters. We compute the 
effective observable. It depends linearly on the improvement parameters. To order 
s + 1 it has no other choice. Then we claim invariance to obtain a linear system of 
equations for the improvement coefficients. This system turns out always to have a 
unique solution: the improved observable. 
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4 Hierarchical fusion rules 



The outcome of our analysis is a system of observables O n ((f>\f) on the renormalized 
trajectory parametrized by /. Under a hierarchical renormalization group transfor- 
mation O n (4>\f) is multiplied by the moving eigenvalue e n (/?(/)) and the coordinate 
is changed to (3(f). This scheme is obviously iteratable which is the reason why we 
introduced it from the beginning. To compute general correlation functions we need 
one more ingredience which is the notion of hierarchical fusion rules. The general 
form of our observables in terms of \P = \<f>\ — r is 



r=0 



n . 



o^m = o, 

m „ 

m=0 

= E <Ih-2- 2 ™ :^- 2 - 2m : . (46) 

m=0 

The normal ordering covariance is 7' = I7. Let us put 7 = 1 to simplify the 
notation. To zeroth order in / we rediscover normal ordered monomials in Their 
perturbations along the renormalized trajectory prove to have no first order terms 
in the running coupling /. To every order of perturbation theory in / we find only 
finitely many normal ordered powers of \E* '. The highest power is n + r — 2 for r > 2. 
Associated with this system of observables is a system of fusion rules defined by 



oMf)o m (<f>\f) = J2 N n,mAf)OiW), 

1=0 

00 

Nn,m,(f) = E<^r- (47) 
r=0 



Furthermore from the fusion rules we obtain a symmetric bilinear form on the linear 
space of observables. It is defined by 

(OMf),O m (4>\f)) = Nn,rn;0U)- W 

The physical significance of this bilinear form is that in the thermodynamic limit 
only the overlap of an observable with the constant term is expected to survive. 
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Since correlation functions are expected to decrease with distance the spectrum 
consists of eigenvalues strictly smaller than one on the renormalized trajectory. To 
zeroth order of perturbation theory we recapture the fusion rules of normal ordered 
products 



^(0 = rilml v 
nml n'\m'\V\ K ' ' ' 

n' = -(l + m — n), 
m! = ~{n + I — m), 

I' = -(n + m- I), (49) 

for \n — m\ < I < n + m and n + m — I 6 2%, zero else. Furthermore to zeroth order 
the observables are orthogonal with respect to the bilinear form 

<^ = < m m!( 7 ') m (50) 

The simple pattern (f49"D, ([SO]) becomes perturbed as one moves away from the ul- 
traviolet fixed point on the renormalized trajectory. The perturbation expansion for 
Ni t i-i(f) to fifth order in / is for instance 

Ar M;1 (/) = -2(N -l)f+ [ 17{N ~ 1)2 - 26(N - 1)) f + 0(f), 
iV ljl;2 (/) = 1 + 0(f), 

iW/) = 24( ^~ 1) / 5 + Q(/ 7 )- (51) 

The perturbation expansion for Ni t2 -,i(f) to fifth order in / is given by 



Ni, 2 -M) = -4(iV 



1)/3+ JM^_M^I)) f + (f) 



NiMf) = l-^^f 4 + 0(f), 

NiMf) = ^l-Vf + 0(f). (52) 
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All other fusion rules are zero to fifth order. We observe that orthogonality is 
violated to fifth order for the first and second observable. 

With the fusion rules we can compute correlation functions of our observables. 
Let us consider for example a general two point function. It depends on the hier- 
archical distance k and the hierarchical lattice size k! of the system. The explicit 
formula is 

k oo k'—k 

(O n ((f>\f)O m ((f>\f)) = II (W)) e m (p j (f)))j2 N n,m;l (/?*(/)) II « (^V)) 

3=1 1=0 i=l 

(53) 

for the two point function on the point of the renormalized trajectory parametrized 
by /. This formula simply says that each observable is renormalized independently 
fc-times. Each renormalization step produces a factor given by the eigenvalue of 
the corresponding observable at the location on the renormalized trajectory. After 
k steps the observables land in the same block and are there fused together. The 
result of fusion is then renormalized (k' — fc)-times to obtain the value of the two 
point function. The thermodynamic limit corresponds to k' = oo. The formula (|53| ) 
is still true in this limit provided that one changes to a different parametrization of 
the renormalized trajectory at the point where the running coupling / diverges. 



5 Numerical results 

When dealing with perturbation theory it is natural to question its validity. To tackle 
this problem we have determined our perfect observables and their corresponding 
eigenvalues numerically. The main technical task is to compute the transformations 

KVW) = -2m(y d/i 7 (C)e-^ +c A , (54) 

tvKOW = /d/i7(c)e ^ + o ■ ( 55 ) 

Both equations can be reduced to integrals of the type 

I F W)= /d/i 7 (C)F(|^ + CI) (56) 
with certain scalar functions F(tp). A shift £ — > ip + ( in the fluctuation field yields 

Z F (H) = (27r 7 )^ 2 / d^Cexp (-i-[C 2 + + i^c) F(\C\)- (57) 
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By using polar coordinates and integrating out the polar angles we find for N = 3 



M\1>\) =*J ^exp (-i-[i? 2 + |^| 2 ]) sinh (^j (58) 

In this form the integral can be evaluated by standard numerical methods. 

Equipped with integrators for ( 54 ) and ( [55]) the strategy goes a follows: The first 



step is to determine Vrt- For this purpose we start with the perturbatively improved 
perfect action Vj^ rt (<^>|r ) at a given radius r as bare potential and iterate the RG 
transformation fl54"D . In each step the resulting potential is driven closer and closer 
towards the RT. After 10 steps we end up with a good approximation of Vrt(0 at 
some radius r = r(r ). In this iterative process the potential V(ip) is represented 
as a cubic spline with N v equidistant knots tpi in the range 1^ = [(p m i n , ^max]- F° r 
each iteration one has to evaluate V = 1ZV at these points. The interval 1^ must 
to be chosen in such a way that for the computation of V'(<p G I v ) contributions 



V((p ^ I J) are numerically negligible in (54) 



In the second step the eigenvalues and eigenoperators of the linearized RG trans- 
formation C^rTZ at this very potential Vr,t(0 are computed. The operator £rt Tr- 
acts on the infinite dimensional space of observables. Naturally a computer can only 
handle the restriction of CrtTZ to a finite subspace. Our program uses the space 
spanned by the operators B m ((p\r) = (|0| — r) m with m — 0, . . . , M. To obtain the 
representation matrix L n>m of CrtTZ in this basis the image of every B m ((j)\r) under 
£rt7£ is numerically expanded by a finite difference method in terms of B n (<f)\r') 

M 

(£ RT n)B m (<P\r) = E L n;m B n {cf>\r'). (59) 

n=0 
M 

Finally the desired eigenvectors O m (r') = Qn,m( r ')B n (r') and their eigenvalues 

n=0 

e n (r') are calculated from L n>m . Fig. |l| shows the flow of the largest eigenvalues 
e m (r'). As in the case of the potential the perturbative predictions are in excellent 
agreement with the numerical results down to a radius of about r' ~ 4. Then non- 
perturbative effects show up forcing the eigenvalues to become smaller. A similar 
behaviour can be found for the expansions coefficients Qn,m{ r ') except that devia- 
tions already show up at r' w 10. 

At r = r cr w 2.04 the effective radius r'(r) vanishes and r ceases to be an 
appropriate parametrization of the RT but the eigenvalues and eigenvectors continue 
to flow against their fixpoint values at the HT-fixpoint. 
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A Hierarchical renormalization group 



In a hierarchical lattice field theory in the form advocated by Gallavotti |G| and by 
Gawedzki and Kupiainen [|GK]1 the massless free covariance (— A) -1 is replaced by 



the hierarchical approximation 



n 



v(x,y) = 7 ]T L m ^5 [L - mxUL - my] . (60) 

m=0 

Its parameters are a positive real number 7 and the positive integer valued block 
size L. The bracket [L _m x] denotes the integer part of L~ m x. The lattice is 
{0, 1, . . . , L n — 1} D with n G IN U {00}. The main property of (|B0"D is the hier- 
archical splitting 

v = L 2 - D C T vC + T, (61) 

C being the unnormalized block average C{x, y) = S x \L-i y ] and T the diagonal part 
T(x,y) = jS X) y. The hierarchical covariance on the right hand side is understood 
to live on a lattice shrunk by the blocking scale L. The splitting (|H]) induces a 
decomposition of the lattice field 4> into a fluctuation field £ and a background field 
ifi. The associated block spin transformation is 

™>W = AT [ dMC)e~ v(Ll ~ ?C7 ^ +c) . (62) 



e 

Here d/xr(C) is the Gaussian measure on field space with mean zero and covariance 
r. The n-fold iteration of fl62|) computes the hierarchical partition function (when 
the normalization is M = 1). The transformation ( |5"2"[ ) preserves locality. For V(<p) = 
Y^x V{4>{ x )) we find 7^(0) — J2x 'RV{(f){x)) in a selfexplanatory notation. (JZ is not 
meant to be linear.) (|62T) then reduces to the non-linear transformation 



e -nvw =Af y d^Qe-viL 1 -^^ . ( 63 ) 

on functions of a single variable <f>. Note that C T il) is constant on blocks and that 
L D is the number of sites in one block. The linear block spin transformation for 
observables is 



/d/xr(C)e- v ( Ll ^ cT ^) 

The expectation value of an effective observable in the effective theory equals that 
in the original theory. (K5T) preserves the property of factorization provided that the 



16 



interaction is local. A local observable 0{<p(x)) is defined as a function of the field 
at a single site x. Local observables are again transformed into local observables 
according to 

= (65) 

The transformation ([54]) is the linearization of (|62| ) at the interaction V((f>). The 
local transformation ( |65D is the linearization of (|63D at the local interaction V{(f) 
divided by the block volume L D . The expectation value of a local observable is the 
72,-fold iteration of (|65|) evaluated at zero external field. The thermodynamic limit 
n = oo follows from an infinite iteration. 

The computation of higher correlation functions requires another ingredience 
called fusion. Let us consider for instance a bilocal observable of the product form 
Oi((f>(xi))02(4>(x2)) depending on the field (f> at two sites x\ and x-i- The hier- 
archical distance of x\ and x% is given by dist(xi,x 2 ) = L k ( Xl > X2 \ where k(£i,x 2 ) 
is the smallest positive integer such that [L~ k( - Xl ' X2 ' > xi} = [L~ k ( Xl ' X2 'xi\. Both ob- 
servable are then transformed independently k(xi, x-^) -times according to the local 
transformation fl5"5] ) until x\ and X2 land in the same block, whereupon the effective 
observables are multiplied together. We call this process fusion. The result is then 
again transformed n — k(xi, a^-times according to the local rule to give the 
hierarchical two point function. In the thermodynamic limit n = 00 the observable 
has to be renormalized an infinite number of steps after fusion. The process can be 
visualized fork. 

Both ( |63|) and (pq ) make sense when the dimension parameter D and the scale 
parameter L are continued to real values. 
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Figure 1: Flow of the eigenvalues e m (r') and the expansion coefficients 
Qn,m{f') along the renormalized trajectory. From top to down ei, e 2 , £3 
and <3o,i> <3o,3) Qo,i- Solid lines represent numerical data, dashed lines 
perturbative results. Calculations were done with N v = 200, y? min = 
0,c/w = 30,M = 10. 
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